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We present a small-but-sizeable magnetic polaron picture where transport at high temperatures 
is activated while at low temperatures it is band-like. We show that both double exchange and finite 
bandwidth effects are important to understand colossal magnetoresistance as well as the coincidence 
of the metal-insulator and the ferromagnetic transitions in manganites. The magnetic transition is 
explained using band-like motion of the polarons. 

PACS numbers: 71.30.+h, 71.38.H-i, 72.20.My, 75.10.-b 



Studies on pcrovskite manganites of the form 
Ai_ 5 B 5 Mn0 3 (A=La, Pr, Nd, etc.; B=Sr, Ca, Ba, etc.) 
have yielded a variety of rich phenomena as a func- 
tion of doping ^-coexisting orbital ordering and layered 
antifcrromagnetism at low doping, simultaneous metal- 
insulator (MI) and paramagnetic-ferromagnetic transi- 
tions at intermediate doping (8 ~ 0.2 — 0.4), and charge 
ordering at higher doping (6 ~ 0.5) Q]. To explain the 
magnetic ordering, de Gennes Q some time ago had pro- 
posed double exchange mechanism wherein, on account 
of strong Hund's coupling between the spin of a mobile 
hole and the spin of the localized electrons, the hopping 
integral of the itinerant hole is reduced by half of the co- 
sine of the angle between the 3/2 spins of the localized 
electrons on neighboring sites. However it was recog- 
nized by Millis et al. || that this mechanism itself is not 
sufficient to explain colossal magnetoresistance (CMR). 
Millis and co-workers have proposed a model Q which 
uses Jahn- Teller coupling between electrons and nuclei. 
However this model does not seem to yield satisfactory 
results away from half-filling. Roder et al. have also em- 
phasized the importance of Jahn- Teller coupling in un- 
derstanding these manganites §. Earlier on Dc Teresa et 
al. j(| have reported evidence for sizeable magnetic po- 
larons above the ferromagnetic transition temperature. 
Recently worledge et al. Q have demonstrated that their 
high temperature resistivity data fits well to an adiabatic 
small polaron model. All in all there is growing evidence 
for a small magnetic polaron picture to explain CMR. 

In this paper we study CMR phenomena in perovskite 
manganites (for i5 ~ 0.2 — 0.4) by considering the car- 
riers as adiabatic small-but-sizeable magnetic polarons 
whose high temperature behavior is hopping type and 
low temperature behavior is metal-like. Our notion of 
small-but-sizeable polaron is similar to that of a nearly- 
small polaron introduced by Eagles Our treatment 
builds on the work of Gosar ||, who included finite 
band width effects by actually considering a small-but- 
sizeable polaron whose wavefunction extends to the near- 
est neighbors. Due to Hund's coupling the electrons on 
the nearest neighbors also have the same polarization. 
Our model includes effects due to electron-phonon cou- 



pling, on-site Hund's coupling between itinerant holes 
and localized electrons, coupling between nearest neigh- 
bor localized spins, and strong on-site repulsion between 
two itinerant holes. To understand MI transition we sim- 
plify the Hamiltonian by accounting for the Hund's cou- 
pling through the double exchange hopping term (|] . We 
show that in the absence of a magnetic field the double 
exchange and finite bandwidth effects can make the MI 
transition coincide with the magnetic transition. We also 
find that in the presence of a magnetic field both double 
exchange and finite band width effects lower the resistiv- 
ity and shift its peak to higher temperatures and thus 
can lead to CMR. 

Within a mean-field approach, we calculate the magne- 
tization (M) of the localized spins as well as the polariza- 
tion (An) of the itinerant hole spins. The magnetization 
M is a function of the external magnetic field and, due 
to the strong Hund's coupling, of also the polarization 
An. The polarization An is in turn a function of Al and 
the external field. Thus the two (M and An) are cou- 
pled and have to be solved simultaneously for a given 
density of holes. We have studied the magnetization at 
different doping values, for both with and without ex- 
ternal magnetic fields, and found that our M values are 
qualitatively in agreement with experimental results [fl0| . 
Furthermore our magnetoresistance values also compare 
favorably with experimental ones jlO) . 

Our starting total Hamiltonian is given by 
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In the above equations Cj i<T (a$) is the hole (phonon) de- 
struction operator, t is the hopping integral, (ij) corre- 
sponds to nearest neighbors, lu^ is the optical phonon 
frequency (h — 1), M q is the hole-phonon coupling, Jy 
is the strength of the spin coupling between neighbor- 
ing localized (S=3/2) spins, K B gives the Hund's cou- 
pling between localized spins and itinerant hole (er = 1/2) 
spin , U is the strength of the same site repulsion, and 
rij = c'j a Cj^ a . Furthermore the H p h part corresponds to 
assuming a single orbital per site which on account of 
Jahn- Teller splitting may perhaps be justified. 

To study transport we use double exchange modifica- 
tion and take the total Hamiltonian to be 



Ht = t D E 4 C 3 + Yl 



(ij) 



+ ^c] Cj . e ^M^ + at f ), (4) 



where t DE = ty/(l + M 2 /M'£)/2, and M s is the satu- 
rated magnetization. Now the mobility is given by the 
Einstein relation fi = q e D(3 where q e is the electronic 
charge, D the diffusivity, and f3 = l/k B T. Includ- 
ing finite band width corrections, as done variationally 
by Gosar |l, to calculate the hopping-regime diffusivity 
Dhop — a 2 /(6r) we obtain the scattering lifetime t for a 
narrow phononic band to be 



l/(6r) = F [A , NA] exp [-20tanh(/?wb/4)], 



(5) 



where a is the lattice constant, wq is the Debye fre- 



quency, 9 = 7 2 



1 



with z being the coordina- 



tion number and, for N lattice sites, 7^0 = M^N 1/>2 . In 
the present analysis the small-but-sizeable polaron con- 
dition is {z t 1 I t h < 1. In Eq. (|), F A and F NA are the 
prefactors for adiabatic and non-adiabatic cases with 



Ft 



^0 [1 + {M/Msf 
2tt 2 



and 



F 



[NA] 



_ t 2 DE ^nsmh(l3u Q /2) 



CJ 



(6) 



(7) 



where the adiabatic term F<m is obtained using argu- 



ments similar to those put forth by Holstein |11|. The 
crossover from the non-adiabatic case to the adiabatic 
case occurs when F[n A ] > F[ A ]. For experimentally stud- 
ied systems the adiabatic regime is of interest. Further- 
more, it should be noted that we need 20csch(/3wo/2) >> 
1 for Eq. (^|) to be valid (for a justification in the narrow 
band case 9 = r y 2 see Ref. ]H|] ). As for the diffusivity 
for band conduction, it is obtained by extending Gosar's 
work Q and calculating the polaronic band energy ]l2] | 



Ej: — 2toE exp [—9 coth(/3w /2)] cos(fc/a) 

1 



(8) 



The above expression is similar to the result due to Ea- 
gles H . Then the diffusivity for band conduction is given 
by 



D ba nd = {\VEA 2 T) =6ra 2 ^ 



(M/Msf 



(9) 



where t = iexp [— 9coth((3u>o/2)]. Then based on Fried- 
man's work ]l3| we take the total mobility (//t) to be the 
sum of the band mobility and the hopping mobility and 
hence the total resistivity (1/ p — ChqePr) to be given by 



4n 



c h q 2 a 2 p 
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_ 2 0csch(^)] + 
exp[-20tanh(^)] 
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where Ch is the density of holes. Here it should be men- 
tioned that, even if t and Kh are of the same order of 
magnitude, we can have t << K B so that double ex- 
change holds. 

To proceed further one needs to obtain the magnetiza- 
tion as a function of temperature. To this end we consider 
the following thermally averaged Hamiltonian 



H„ 



H(t) + H, 



(11) 



Next we note that | J y | << K H and that U » t {&K H ) 
and hence completely project out double occupation (see 
Ref. for details) . From the above Hamiltonian H mag 
it follows that, within a mean- field treatment, the mag- 
netization in the presence of a magnetic field H is 



,5' 



M_ _ Es z S z exp[-($/2 + gf i B H)S z (3] 
' M S ~ Es, exp[-($/2 + gp B H)S z (3] ' 



(12) 



where $ = K H {n\f(n l ) - n 1 /^)}, f{n a ) = 1/(1- 
n Q ), and n a is the probability of occupation of a site by 
spin a hole and is given by 



n a = -Yn 

N 4^ 



exp[-/3(*CT + i u)] + r 



where * = (K H SM/M S + gp B H)f(n a ), a = l/2(-l/2) 
for spin a =f (J,) holes and 



e2 = 2th(S)f(n a ) 



1 



cos(fcfa) << k B T c 



(13) 



In the above equation h(S) = 1 — 5 and Tq is the ferro- 
magnetic transition temperature whose value, by treating 
M and An = n" 



n- 1 as small parameters in Eq. (|12|), is 



obtained to be 



knT, 



BJ-C 



'20(1 



9(2-5) 



%K H S\o\ 



(14) 
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We see that Tc increases with increasing S for < S < 
2/3 and that it is independent of both t and Jy. Fur- 
thermore because of Eq. ( |l3| ) the values of M, n a , and 
Tc are all independent of dimensions [see Eq. (|l|)] . Here 
it must also be mentioned that when double occupancy 
is allowed h(6) = /(n T ^) = 1. 

Using the constraint that rv + nr = S, we can obtain 
An and M by solving Eq. jl2] ) . In Fig. |TJ we have plotted 
the magnetization ratio M/Ms as a function of the re- 
duced temperature T/Tq for S = 0.3 and 0.4, g = 2, and 
magnetic fields H = 0T and 15T. We have assumed a 
smaller value for the Hund's coupling (Kh ~ 0.08586^) 
than what seems to be its value based on experiments 
(~ XeV) because we wanted to set Tc — 300A" at 5 = 0.3. 
Alternately one can also get lower Tc by assuming, as 
suggested in Ref. ||, that only a small fraction of the 
dopants yield mobile holes. The values of the magneti- 
zation for H = 15T at Tc are sizeable because of the 
tendency of the system towards a ferromagnetic phase. 
Here it should also be mentioned that An attains sat- 
uration values much faster than M fl5| . We have also 
calculated the magnetization curves with double occu- 
pation of a site being allowed and find that the M/Ms 
values for with and without double occupation being al- 
lowed are close to each other both in zero field and at 
15T [ p^[ . Although our magnetization curves are quali- 
tatively similar to the experimental curves of Urushibara 
et al. jTTJ], the experimental M/Ms values rise faster as 
T is lowered. 

We will now discuss the resistivity given by Eq. (10). 



The conduction goes from a hopping type at high tem- 
peratures to a band type at low temperatures. In Fig. ^ 
we have shown the dependence of resistivity p on tem- 
perature at various magnetic fields. The general trend of 
the resistivity including the drop at the MI transition at 
H = 0T is similar to the experimental results Jl(]]. On 
introducing a magnetic field the system gets magnetized 
at temperatures higher than Tc and thus the value of 9 
is smaller (see Eq. jicj)). Consequently the resistivity is 



smaller and T p max (the temperature at which resistivity 
becomes maximum) increases Jl6|| , 

For T > Tc, when D^and/ T>i lop >> 1 the magnetore- 
sistance Ap/p(0) = (p(H) — p(0))/p(0) is given by (see 
Eq. ©) 



Ap/p(0) S3 exp 



(z+1) t 2 M 2 pu^o, 
2 7 2 co 2 M 2 [ 2 ' 



1, (15) 



and when D band /D hop « 1 it is given by 



Ap/p(0) 



exp 



(z + l) t 2 



^tanh(^) 



1 + (M/Ms) 



1. 



(16) 



For a fixed value of the reduced temperature T /Tc, an 
increase in the ratio psH/Kn increases M/Ms and con- 
sequently the magnetoresistance also increases. 



Actually T p ma?> (the temperature at which the resis- 
tivity given by Eq. (10), after taking M = 0, attains a 
maximum) need not be equal to the ferromagnetic tran- 
sition temperature Tc- If Tpma* < Tc, by decreasing 

7 2 or increasing activation energy (8ujq/2) decreases 

and T p max can be increased Jl - 



to be made equal to Tc 
and this also increases the magnetoresistance (see Eqs. 
@, ([TJ), and (jlj)). For T p mo» < T c , the MI transition 

can still occur at Tc if yr is sufficiently large |IJ| 

while if 2y ,j! Z7 1 * s ver y small the MI transition occurs 
below Tpmax as can be seen from Eq. (|l0|). The other 
case, where T p max Q > Tc, corresponds to MI transition 
occurring at a higher temperature than Tc and is in any 
case not experimentally observed p5| . 

In Table | we report the calculated values of magne- 
toresistance —Ap/p(0) at Tc and the optimum values of 
7 2 (obtained when T p maa = Tc) for doping 5 equal to 
0.3 and 0.4, Debye temperature Tq = 500A, and for var- 
ious values of the dimensionless hopping integral t/uiQ. 
We find that the magnetoresistance increases with in- 
creasing values of t/u>o thus showing the importance of 
bandwidth. Also 7 2 pt values increase with increasing t/uo 



because T, 



n,a:r 
1 M = 



Tc- Furthermore, it is mainly due to 



the larger values of (M/Ms) f° r 8 = 0.4 compared to 
those of (5 = 0.3 that the values of — Ap/p(0) are larger 
for 5 = 0.4. It appears that our model can give magne- 
toresistance values comparable to the experimental ones 
JTo| | . In fact one can get a larger magnetoresistance by 
taking a smaller Tc value but keeping ojq/ (ksTc) fixed 

From Eq. @ (or Eqs. @ and @) we see that 
for small values of M /Ms the magnetoresistance (for 
T > T c ) is of the form -Ap/p(0) = C(M/M S ) 2 where C 
is a constant of proportionality. We found, for the cases 
considered in Table 1, that the optimum values of 7 2 
that make 2p™°L x = Tc are such that Dband/Dhop < 1 so 
that the magnetoresistance can be qualitatively given by 
Eq. @. In Eq. @, close to T c , tanh(/3u /4) ~ AV 4 - 
From Eq. (0) we see that Tc increases with the doping 8 



and hence the constant of proportionality C (oc 



(2-<5) 



decreases with increasing S which agrees with the findings 
of Ref. (h]]. Furthermore the coefficient also increases 
with increasing values of it^ . We have calculated 
values of C by treating M/Ms as a small parameter in 
the exact expression for —Ap/p(0) at T — Tc- When 
t/uio = 4(8) and 7 2 = 9.5(15.0), for <5 = 0.3 we get 
C fa 3.8(8.5) while for S = 0.4 we obtain C « 3.1(6.4). 
Our calculated values of C are larger than those reported 
in Ref. |l(J. Past attention [ p~7|JT^ ] has focused at de- 
pendence of C on the ratio Kh /t in Kondo lattice type 
models that ignored electron-phonon coupling. While In- 
oue and Maekawa jlTj for Kh — > oo obtained C = 7/4, 
Furukawa (Til found that the value of C increased with 
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increasing values of Kn/t and that at larger values of 
Kff/t the value of C decreases with increasing doping. 

In conclusion we say that both double exchange and fi- 
nite band-width corrections are important to understand 
CMR. In our picture, adiabatic small-but-sizeable mag- 
netic polarons are involved in activated transport at high 
temperatures and metal-like conduction at low temper- 
atures. At the MI transition, the band-like motion of 
the carriers also produces a paramagnetic-ferromagnetic 
transition due to strong Hund's coupling between itiner- 
ant and localized spins. Studying the transport behavior 
at low temperatures, including a Fermi liquid analysis, 
is left for future. The effect of including both d 3z 2_ r 2 
and d x 2_ y 2 orbitals also needs to be investigated for our 
model. Lastly we note that as the system's temperature 
is lowered below Tc the magnetization increases and con- 
sequently the activation energy (Oluq /2) decreases and the 
polarons tend towards large polaronic behavior. 
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FIG. 1. Plot of the magnetization ratio M/Ms versus 
the reduced temperature T /Tc when no double occupa- 
tion is allowed, doping 8 = 0.3 (and 0.4), Hund's coupling 
Kb « 0.0858eV, and magnetic fields H = 0T and H = 15T. 



FIG. 2. Plot of the resistivity p in units of 4n/(chq 2 a 2 ) 
versus temperature T in 3 dimensions when no double occu- 
pation is allowed, 8 = 0.3, Kh ~ 0.0858eV, dimensionless 
hopping integral t/uio = 6, optimum j 2 = 12.2, Debye tem- 
perature Td = 500A", and for the following magnetic fields: 
(i) H = 0T; (ii) H = 15T; (iii) H = 30T; and (iv) H = 45T. 



TABLE I. Calculated values of the magnetoresistance 
— Ap/p(0) at Tc and the optimum 7 2 for various values of 
t/uJo, To = 500A', 8 = 0.3 and 0.4, magnetic field H = 15T, 
K H w 0.0858eV, and T p ™* = T c . 
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